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Abstract. A dynamical price formation model for financial assets is presented. It 
aims to capture the essence of speculative trading where mispricings of assets are used 
to make profits. It is shown that together with the incorporation of the concept of risk 
aversion of agents the model is able to reproduce several key characteristics of financial 
price series. The approach is contrasted to the conventional view of price formation in 
financial economics. 



INTRODUCTION 

A remarkable feature of (financial) economics is that an immensely complex 
system - consisting of a relatively large number of human beings, interacting 
through the exchange of goods, information, services, and money - reduces to 
one-dimensional functions in time, the prices of assets. These prices potentially 
carry substantial information about the underlying system, and are therefore an 
intriguing subject to study. 

Technically the asset price, be it of a stock, option, or any other liquid financial 
derivative can change whenever a trade (interaction of market participants) takes 
place. Prices are formed through the actions of three types of agents: A seller 
offers an asset at a certain price, usually through a sell order to the market maker. 
Potential buyers offer to buy the asset at a certain price, or at any price. The market 
maker sets the actual unique price at which trades then occur. For convenience 
one frequently works with the log-return process rather than with the price process 
itself. The log-return over a specified time interval (aggregation time r) is simply 
given by the logarithm of the quotient of prices separated by time r: 



r.(t) = log(p(t + r))-log(p(t)) . (1) 

Due to the addition property of logarithms it is easy to change from a given aggre- 
gation time r to higher ones by taking the prices which form the boundary of the 
time interval of interest. These return series possess a highly non-trivial structure. 



and over the last decades considerable effort has been devoted to quantitatively 
characterize such processes. 

Maybe the best known fact about financial price time-series is that the distri- 
bution of returns is leptocurtic. On short time-scales (r ranging from minutes to 
several days), for small price fluctuations the return distribution can be well ap- 
proximated by a Levy stable distribution [1] , while its tails - emerging from violent 
price changes - generally follow a power law, see e.g. [2]. It is well established 
that for longer time-scales, i.e., aggregation times r of three weeks and longer, the 
distribution functions of returns are almost Gaussian. Historically there has been 
some confusion on this subject: In the 1960s Mandelbrot and Fama found empir- 
ical evidence that the cumulative distribution functions of returns of time-scale r, 
/(r^), behave like stable Levy distributions, which are characterized by a parameter 
1 < /i < 2, where fi = 2 corresponds to the special case of a normal distribution. 
For 11 < 2 the distributions are fat tailed and higher momenta - including the sec- 
ond - are infinite. They estimated an empirical value oi ^l — 1.7, indicating that 
higher momenta will not exist and therefore putting in doubt most statistical quan- 
tities. In the 1980s power laws for the tails of the distributions have been found, 
i.e. f{Tr) ~ krl""; with a > 2. This result has been confirmed throughout the 
1990s and the exponent is nowadays believed to be close to a = 3 [3]. 

For hquid assets, C(100 — 1000) trades per day, the stock returns are basically 
uncorrelated, i.e., the autocorrelation functions with time lags larger than one are 
usually compatible with zero. However, one notices the fact that there are periods in 
time when price fiuctuations are stronger than at other times. This is called volatil- 
ity clustering and can be captured quantitatively by the autocorrelation function 
of the absolute values (or squares) of the returns. These volatihty autocorrelations 
do not vanish exponentially but again show an empirical power-like decay, lag"''', 
with /3 ~ 0.1 - 0.5, see e.g. [4-6]. 

There have been observations that higher momenta of volatility (Ir^l"^) are scaling 
processes with respect to aggregation time r and that the corresponding scaling 
exponents are not linear functions in g, which is an indication for a multiscaling 
process [4,7]. 

Further empirical features of price processes are that prices are non-stationary, 
the second moment of returns (r^) exists (but converges slowly), whereas the fourth 
moment {r*) does not. Prices are usually quantized (1/16 dollar at the NYSE), 
and can not jump arbitrarily much from one trade to the other. The closing price 
is the price of the asset at the end of a trading day. 

In this paper we model the one-dimensional price process, by casting some fun- 
damental characteristics of economic systems and human behavior into a set of 
differential equations, and compare the emerging price time-series statistics with 
known facts of real price series. The paper is organized in the following way: The 
next section is a mini-overview on the standard views in (micro) economics of price 
formation. We then motivate and propose a dynamical system, show some of its 
features, and discuss the results. 



STANDARD TREATMENT IN ECONOMICS 



Financial economics literature, in particular the field of micro-economics, is dom- 
inated by several concepts which are powerful from a theoretical standpoint but 
which are known to have their limits when it comes to the analysis of real time- 
series. These concepts are: 

Dogma 1: Prices are the equilibrium result of trades, i.e., all sellers and buy- 
ers come to an optimum state (Parcto optimum), where everybody is better off 
than without doing the trades. It can be shown quite generally with the help of 
the Brower fixed point theorem, that such an equilibrium state exists which also 
provides an optimum [8]. For this to be true it is important that the influence of 
trading barriers is kept small. In this approach all the agents involved in trading 
maximize their von Neumann- Morgenstern utility functions. These are monoton- 
ically rising convex functions, which describe the "happiness" of the individual 
agent as a function of his wealth. The fact that utility functions rise monotoni- 
cally reflects that "more is preferred to less" , the convexity takes into account the 
risk-aversion of agents. The main problem associated with this view is that it is 
not able to account for speculative trading at all. On the other hand one knows 
that speculative trading is dominating entire markets. For example in 1995 the 
daily trading volume at the foreign exchange markets exceeded lO-*^^ USD which 
was about 20 times the daily world gross national product. 

Dogma 2: Markets are efficient. The basic requirement for efficiency is that all 
participating agents are fully rational, that there is rapid information flow, there 
are no transaction costs or other market friction, and that everybody has identical 
access to all the relevant information. As a consequence, prices simply reflect the 
expectations on future earnings from investments made today. Thus trading would 
only occur at the appearance of news which would influence those expectations. 
Other than that there would be no trading activity and speculative trading would 
be entirely useless. Another consequence of the efficiency hypothesis is that there 
exist no patterns in financial time-series which could be exploited in one or the 
other fashion. Also there should not exist arbitrage (risk-less profit) opportunities. 
In fact the non-existence of arbitrage forms the basis of many celebrated results in 
financial mathematics, especially in the field of pricing financial derivatives. 

Obviously the efficient market hypothesis is rather unrealistic and there are 
strong problems related with it: Agents are not fully rational, but often tend to 
have non-rational strategies, such as heard behavior, believe in experts' opinions, 
etc. Moreover it has been clearly demonstrated that most trading is not infiuenced 
by news [9]. Finally, it is needless to mention that there exist financial arbitrage 
companies which are specialized in exploiting patterns in financial time-series. 

Dogma 3: Mathematical tools: In financial literature the most commonly used 
mathematical tool is game theory. Most of the time, non-iterative one or two 
period games are considered, which severely limits any study of the dynamics of 
a system. If dynamics is studied, usually the corresponding models are based on 
Wiener processes (Brownian motion). For example as a starting point for many 



models a price process is assumed to have the form 



dp{t) = II pit) dt + a pit) dWit) 



(2) 



with jji the so-called drift term, a the volatility, and Wit) the Wiener process. The 
first term on the right side determines the overall exponential growth of an asset 
price with growth-rate /i. Superimposed on that, the second term introduces a 
source of stochasticity, whose relative importance is controlled by the volatility a. 

Wiener processes are nowadays mathematically well understood and are rela- 
tively easy to handle, since they are Gaussian processes. These kind of processes 
have led to a number of celebrated results in option pricing [10], interest rate mod- 
els, etc. For an overview see [11]. However, there are severe problems associated 
with Wiener processes if one tries to explain real price series. Clearly, in Gaussian 
price processes no power laws are present in the return distributions, there is no 
volatility clustering and there exist no power laws in the volatility autocorrelations. 

In order to make things more realistic, so-called ARCH-GARCH processes have 
been introduced [12]. Such models are nowadays frequently used in everyday eco- 
nometrics and risk management. A G ARCH (1,1) process has one source of ran- 
domness and obeys the following evolution equations: 



where Zt arc uncorrclatcd normal random variables, and a, h and c arc real numbers. 
Such models are able to explain clustered volatility, but fail to give power laws in 
the volatihty autocorrelations. The main problem with GARCH models is that 
they are ad hoc models and do not relate the variables to an economic context. 

Obviously there is a need for a better understanding of the origin of the basic 
features of price processes. Recently there have been put forward some new ideas 
from the mathematics and physics community, some of which are certainly more 
successful in describing reality than the standard methods used in economics. Just 
to mention a few directions of new developments, there are contributions from 
evolutionary game theory [13], agent based models and minority games [14], and 
spinglass-inspired models [15]. Stochastic models of market maker behavior have 
been discussed in [16] and in market impact models [17] several realistic price for- 
mation scenarios are studied. In the following we present a dynamical systems 
approach, which tries to model fundamental behavior of speculating market partic- 
ipants. The aim is to capture this behavior in a set of coupled differential equations. 
The solution to this system yields the price process. 




atzt ; Zt ~ M{Q, 1) 



(3) 



THE DYNAMICAL THERMOSTAT MODEL 



The model is supposed to capture three fairly general features of price dynamics 
and investor behavior: 



Fact 1: If money is invested in a bank or in government bonds, which is often 
referred to as a risk-less investment, wealth will increase exponentially with a fixed 

interest rate (for short-term investments). If one decides to invest in risky assets 
like stock, this exponential tendency (drift) should also be present with a some- 
what larger expected rate of return, which compensates for taking the risk (risk 
premium) . The price dynamics of a risky asset will be more complicated since it is 
coupled to a complex dynamical system - the market. This can be written as 

p{t)=np{t)-cUt)p{t) , (4) 

with n the drift coefficient and c a coupling constant, is a coupling or friction 
variable, whose dynamics is the main subject of this work. Note that this equation 
is formally similar to the standard price dynamics using Wiener processes, Eq. (2). 

Fact 2: The essence of speculative trading is that agents try to make profits by 
using mispricings of assets. Suppose there exists the "true" value of an asset at 
a given time %){t). The actual price p{t) can, and in general will, differ from 
this value. As an example of how such a mispricing can be used to make profits 
imagine that the asset is "under-priced", i.e., the actual price is lower than the true 
value, p < 7^(i). Alert agents who notice this difference will buy the asset and hold 
it until the market as a whole (not that alert) values the asset correctly. In this 
course the price of the asset will rise by an increased demand, and the mispricing 
gets reduced. Obviously the rise in price can be used by the agents to realize 
profits. If an asset is "over-priced" there are also ways to make profits by so-called 
"short-selling" . In this context it seems reasonable to view the price changes in 
Eq. (4) proportional to the mispricing: 

i^^\\p-%{t)] , (5) 

with Ti being a constant, which can be interpreted as the on-set time of the cor- 
rectional movements. As long as the asset is under-priced, is negative, will 
decline and eventually fall below zero, which leads to a positive second term on the 
right side of Eq. (4), which makes the price rise. 



TABLE 1. Meaning of model constants and variables and its relation to physics 



symbol 


meaning in model 


relates to in physics 
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price function 


momentum 
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risk-free rate (bank rate) 


"force" term 
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primary coupling of the asset to the market 
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FIGURE 1. (a) Mean (boxes) and variance (circles) of as a function of Ti{t) = 1, 2, • • • , 10. 

In (b) several time-series of Ci(i) are shown with the corresponding histograms, which indicate 
fat tailed distributions. 



Fact 3: Price changes are not un-restricted, which can be explained by the fact 
that agents are risk averse. Agents are aware that they can misjudge the situation 
about mispricings and as a consequence make losses. Since they are not entirely sure 
if their investment will be profitable or not, they will invest only limited amounts. 
The model should thus have a handle to restrict the price changes to a region 
around zero with controllable variance. Technically this can be done by introducing 
a new dynamical variable ^2 and by adding a term —^1^2 to the left hand side of 
Eq. (5). The dynamics of ^2, which keeps the variance of ^1 at a value of 7i, is 
given by 



i2 = ^[ei-m) 
^2 ^ 



(6) 



where T2 is another time constant. Ti{t) is a factor (maybe time dependent) that 
can be interpreted as a measure of collective inverse risk aversion of the involved 
agents. If Ti is small, agents are not willing to take risk and only small price 
changes will occur. By gathering the above arguments we propose the following 
model: 



p = 








ii = 






6 = 


1 







(7) 



The variables are collected in Table 1. The first two equations have the form of 
dynamical thermostat equations which have been introduced some time ago [18]. 
To build the bridge to statistical physics, if one considers instead of p in Eqs. 



TABLE 2. Parameter dependence of log-return dis- 
tribution functions, a corresponds to the tail expo- 
nent at a given aggregation time t. For all simula- 
tions we took 1000 time steps, and fixed 71 = Tq. 
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2.62 
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1.60 


0.30 
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0.0037 


0.120 


1.36 



(7), the first equation stays the same up to a rescaling factor 2, and the second 
equation can be seen as a dynamical thermostat, which keeps the kinetic energy 
(p^) at a temperature 7o(t). ti is then the thermostat on-set time. 

Model features 

The model is written in continuous variables. However, trading is not a con- 
tinuous process but occurs at specific points in time. To capture this feature we 
solve the equations numerically with a Runge Kutta 4th order solver with time 
increments At of 1/1000 - 1/10000. Every time increment is supposed to model 
one single trade (tic-scale). The prices are recorded at integer times 1,2,...,N which 
represent the "closing prices", i.e., one time unit is 1 day. At can thus be consid- 
ered an additional (hidden) parameter of the model. In the following we decided 
to use a dynamical step size, according to predetermined error tolerances. The 
variable of interest is p, which will be analyzed in the following in the identical 
manner as realistic stock data. In Fig. la we show by numerical simulation that 
the third equation in Eqs. (7) controls the variance of price changes, such that 
7^=Var(^i). The mean values of ^i(t) (boxes) are clearly compatible with zero, the 
variance (circles) rises linearly with the (constant) value of 7^. The corresponding 
distributions of the ^i{t) appear to be fat tailed but have the desired variance. The 
time-series and corresponding histograms of friction variables in a generic case 
are seen in Fig. lb. The solid hue in the histograms is exp[— (^i/71)^]. 

For simplicity we set /x = for the remainder of the paper. To model a realistic 
situation, we chose '7o(t) = exp(rt) with r = 0.0007 which corresponds to an 
average growth of about 19% per year (~ 250 trading days). We assume that the 
risk aversion of agents becomes less as their wealth increases (convexity of utihty 
functions), which we took into account by simply setting 7^(t) = Certainly 
one could think of other choices which capture the same feature. The on-set times 
Ti and T2 have been chosen to be less than one, which states that decisions on 
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FIGURE 2. (a) Model price series for 300 time steps, following a To{t) = exp(rt). (b) The 
corresponding return process exhibits volatility clustering which leads to non-trivial volatility 
autocorrelation functions which are compatible with a power decay, (c). In (d) a part of the 
trajectory of the same process is shown. 



mispricings are made on a time-scale below one day, which is nowadays a realistic 
assumption for liquid assets (intra day trading). With the coupling constant c, it is 
possible to weigh the fluctuations relatively to the drift term. In the following we 
set it around 0.1 and 0.3. In Fig. 2 we show model results for the price series (a) 
and the corresponding log- returns (b). Clearly, the price follows the exponential 
form of %{t) (dashed line) and the returns show clusters of enhanced volatility. 

To become more quantitative we looked at the first 100 lags of the normalized 
autocorrelation functions of squared log-returns. In Fig. 2c the situation is shown 
for the same parameters as before. The first 100 lags of the autocorrelation is seen 
to be compatible with a power decay (straight line) with an exponent of the order 
of /3 ~ 0.1, which is within the realistic range. 

In Table 2 we gather the power exponents a of the return distribution functions 
for an aggregation time of r = 1 "day" . These values have been obtained by linearly 
fitting the tails of the distribution functions in a log- log plot, and by averaging over 
3 independent runs. In most cases a lies in the range between 2 and 3, depending 
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FIGURE 3. Comparison of higher momenta of real (left) and model data (right). Top: The 
higher momenta of the return process are scaling processes with respect to aggregation time r. 
Bottom: {\tt\'^) / {\rT-\)'^ is not a constant, indicating a multifractal underlying return process. 

somewhat on the on-set times ri and T2. This result is again in good agreement 
with observed data. 

Finally we computed the higher momenta dr,-!'') of the return processes for vari- 
ous q values, Fig. 3 (top). We observe that (|r^|^) scales with respect to aggregation 
time r. The corresponding log-log slopes do not depend on q linearly, which is in- 
dicated by a slightly declining {\'rT\'^) / {\r-r\Y, Fig. 3 (bottom). This suggests that 
the underlying process is a multifractal. When comparing to realistic stock data of 
about the same data length (left), again, nice agreement is found. 

The presented model has a rich dynamical structure which will be subject for a 
more detailed investigation in future work. Here we just mention that formally the 
stationary solutions are given by p = 0, ,^1 = ±-\/7i, and ^2 = The 
determinant of the Jacobian matrix is Det(M) = — 2c^f /r| which for the stationary 
solutions reduces to — 2c7^^''^/r|. For the same set of parameters as before the 
attractor is plotted in Fig. 2d for several time steps. Note that the attractor is 
on the "tic-scale" and that the price process used for the above analysis takes only 
every 1000*'^ point (on average) from that trajectory. 

DISCUSSION AND CONCLUSION 

We have presented a non-equilibrium price evolution model, which captures reg- 
ulatory price movements through the agents' desire to use mispricings to make 
profits, and which takes into account risk-aversion of those agents. We contrasted 
this model to the standard believes of financial economics where prices emerge from 
equilibrium, by agents maximizing their utility functions. In our approach we do 
not use the standard mathematical tools like Wiener processes or game theory. 



The apparent differences of our dynamical system with the standard view is the 
exphcit use of a non-equihbrium concept, which is reahzed by the introduction of 
thermostats which model the behavior of agents. Our model does not contain any 
sources of randomness, the erratic behavior originates from the non-linear nature of 
the model. Similarities to the standard approach are that we also use the concept 
of risk-aversion and that the basic price equation is formally similar to the standard 
formulation. The presented model - even though very simple conceptually - leads to 
realistic looking price dynamics over a wide range of parameter settings. Volatility 
clustering, fat-tailed distributions of returns, correct looking autocorrelation func- 
tions, higher momenta and multifractal spectra are well reproduced. Of course, a 
purely deterministic model can never be the full truth of a complex system like 
financial markets, but it is intriguing that it is possible to explain most stylized 
facts of its resulting time series. Whether there are practical sides to the model 
remains to be seen. To mention an apparent consequence, in the present framework 
the pricing of financial derivatives would be fundamentally different from current 
practice where arbitrage free prices are derived from replicating portfolios, which 
are chosen such that the (Wiener) random sources cancel each other. 
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